Nuclear Charge Radii of 7,910 Be and the one-neutron halo nucleus n Be 
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Nuclear charge radii of 7 ' 9 ' 10 > 11 Be have been determined by high-precision laser spectroscopy. On- 
line measurements were performed with collinear laser spectroscopy in the 2s%/2 —> 2pi/2 transition 
on a beam of Be + ions. Collinear and anticollinear laser beams were used simultaneously and the 
absolute frequency determination using a frequency comb yielded an accuracy in the isotope-shift 
measurements of about 1 MHz. Combination with accurate calculations of the mass-dependent 
isotope shifts yield nuclear charge radii. The charge radius decreases from 7 Be to 10 Be and then 
increases for the halo nucleus n Be. When comparing our results with predictions of ab initio nuclear 
structure calculations we find good agreement. Additionally, the nuclear magnetic moment of 7 Be 
was determined to be —1.3995(5) (is and that of Be from a previous /3-NMR measurement was 
confirmed. 

PACS numbers: 32.10.Fn, 21.10.Ft, 27.20.+n, 42.62.Fi, 31.15.ac 
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The discovery of halo nuclei in 1985 [l[ triggered a large 
number of experiments in order to understand the prop- 
erties and behavior of these extraordinary systems. In a 
halo nucleus, individual nucleons - in most cases neutrons 
- can reside far away from the nuclear core. The best es- 
tablished and investigated neutron-halo nuclei are 6,8 He, 
11 Li, and u Be. In a neutron-halo nucleus, the change of 
the nuclear charge distribution provides information on 
the interactions between the different subsystems of the 
strongly clustered nucleus. Such changes may originate 
from two effects: firstly, the motion of the nuclear core 
relative to the center of mass, and secondly core polariza- 
tion induced by the interaction between the halo neutrons 
and the core. In the case of helium the nuclear core is 
tightly bound and only a 4% contribution of core exci- 
tation is expected Q. The dominant part of the change 
is induced by the center-of-mass motion. Conversely, the 
9 Li core in n Li is rather soft and core polarization can 
play a much more important role Q • Separation between 
these effects is not trivial since the correlation between 
the halo neutrons needs to be considered. Thus, to de- 
velop a consistent picture of neutron-rich light nuclei, 
more data are necessary for different systems, which is 
a motivation to determine accurately the charge radii of 
7,9,io,iig e ii j$ e i s the first one-neutron halo nucleus for 
which the charge radius is reported. 

The only known model-independent way to determine 
charge radii of unstable isotopes is the laser spectroscopic 



measurement of the nuclear-volume dependent isotope 
shift. Nuclear-volume shifts for isotopes of the light el- 
ements are very small, and short-lived isotopes can be 
produced only in marginal quantities. Therefore, high- 
precision measurements are extremely difficult. More- 
over, it is a major challenge to separate these small effects 
from the mass shifts which dominate by three orders of 
magnitude. A breakthrough on the theoretical side was 
the sufficiently accurate calculation of mass shifts for sys- 
tems with up to three-electrons 0, @, 0, H] • These re- 
sults have been used for the first time to determine the 
charge radii of 6 ' 7 ' 8,9 Li and 6 He (lp| . Later, the charge 
radii of the two-neutron halo nucleus n Li ll| and re- 
cently that of the four-neutron halo nucleus 8 He [l2| were 
determined based on these calculations. The results are 
now benchmark tests for nuclear structure theory. 

Ab initio descriptions of the structure of light nuclei 
based on realistic two- and three-nucleon interactions be- 
tween individual nucleons have been developed during 
the past years. Such are Greens-Function Monte-Carlo 
(GFMC), Fermionic Molecular Dynamics (FMD), and 
Large Basis No-Core Shell Model (NCSM) calculations, 
which have predicted nuclear charge radii of beryllium 
isotopes. Moreover, the charge radius of 7 Be gives con- 
straints for the determination of the 7 Be(p, 7) 8 B iS-factor 
that is crucial for solar neutrino experiments HI, 3] • For 
these reasons, the isotope shift of beryllium isotopes has 
recently attracted considerable interest. Measurements 
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FIG. 1: Setup for collinear laser spectroscopy with parallel 
and antiparallel excitation and a frequency comb for reference. 



on trapped and cooled ions in a Paul trap have been pro- 
posed for the Dl line [l~5[ and first results in the D2 line 
were reported [l6j |. However, these measurements have 
not yet reached the accuracy required to extract the field 
shift. 

We have developed a technique to determine the opti- 
cal transition frequencies in collinear laser spectroscopy 
and used it for on-line measurements on 7 > 9 > 10 > 1:L Be. The 
radioactive isotopes were produced at ISOLDE (CERN) 
with a 1.4 GeV proton beam impinging on a uranium 
carbide target. Beryllium atoms were ionized with the 
laser ion source (RILIS). After extraction, acceleration 
to 50 kV and mass separation, the ions were delivered 
to the collinear laser spectroscopy setup [l8|. This has 
been used for isotope-shift measurements on many ele- 
ments down to neon (Z = 10) [19(. However, for very 
light elements the systematic uncertainties caused by the 
measurement of the acceleration voltage exceed the nu- 
clear volume effect. Now, we can overcome these limi- 
tations employing a frequency comb and measuring the 
absolute transition frequencies v v for parallel and v a for 
antiparallel geometry of ion and laser beams. This yields 
the rest frame frequency vq independently of the acceler- 
ation voltage via the relativistic formula = v v ■ v a . 

A schematic layout of the setup and the laser system 
is shown in Fig. [TJ The ion beam is superimposed with 
the co- and counter-propagating laser beams by an elec- 
trostatic 10° deflector. The laser beams enter and leave 
the beamline through a pair of Brewster windows. Two 
adjustable irises and further sets of fixed size apertures 
down to 5 mm diameter are used along the beam line to 
ensure good overlap and to avoid stray light in the optical 
detection region. Doppler tuning is performed by chang- 
ing the potential of the detection region in the range of 
±10 kV. Two photomultipliers are used for resonance flu- 
orescence detection. 

The output of two dye lasers was frequency-doubled to 



produce ultraviolet light at 313 nm and the UV beams 
were well collimated over a distance of about 8 m with a 
beam diameter of about 3-4 mm, well adapted to the ion 
beam size. To avoid strong optical pumping and satura- 
tion broadening of the induced transitions, the UV light 
was attenuated to less than 5 mW. The fundamental light 
at 626 nm for collinear and anti-collinear excitation was 
frequency-stabilized in different ways. One of the dye 
lasers was locked to an iodine line using frequency modu- 
lation saturation spectroscopy, while the second one was 
locked to a frequency comb (Menlo Systems FC1500). 
During the beamtime, the frequency of the iodine-locked 
laser was repeatedly measured with the frequency comb. 
In total, 12 different iodine lines were used and their fre- 
quencies measured with standard deviations on the order 
of 20 kHz. The Rb clock that was used as a reference os- 
cillator for the frequency comb introduces an additional 
systematic uncertainty of about 350 kHz. However, this 
contribution cancels out for the isotope shifts in which 
frequency differences are evaluated. 

Measurements were performed with the frequency of 
the collinear laser locked to a suitable iodine line. The 
voltage of the detection region was tuned to record the 
hyperfine-structure pattern of the collinear excitation. 
Then, the frequency of the anti-collinear dye laser was 
locked to the frequency comb and its frequency chosen 
such that the resonance pattern was covered by the same 
voltage range. Repetitive scanning was performed by 
tuning the voltage across this range. Dwell times of 20 ms 
per step were used, resulting in about 4 s per scan. Two 
remote-controlled mirrors were used to block alternately 
one of the two laser beams during each scan. Typically 
2 x 50 scans were accumulated for isotopes with high 
yields, whereas up to 2 x 400 scans were taken for a spec- 
trum of 11 Be (see Fig. [2]). 

The spectra were fitted with Voigt profiles of common 
widths for all hyperfine structure components. The hy- 
pcrfine pattern was reproduced by direct calculation of 
the F — > F' resonance position from 
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~[A 2p C(F,I,J)-A 2s C(F',I,J)} (1) 



with the center-of-gravity frequency v cg , the A factors 
of the ground (A 2s ) and excited states (A2 P ), respec- 
tively, and C(F, I, J) = F (F + 1) -I (I + 1)- J(J+ 1). 
Chi-square-minimization was done by varying A and v cg 
amongst other parameters. Weak satellite peaks aris- 
ing at higher acceleration voltages due to inelastic col- 
lisions with residual gases were also included in the fit- 
ting function. The achievable accuracy was tested in an 
off-line beamtime and systematic shifts of the transition 
frequency evaluated. Only a small effect caused by beam 
misalignments could be observed. 

The obtained isotope shifts Su^g — ^o( A Be) — fo( 9 Be) 
are listed in Tab. [TJ The quoted uncertainty represents 
the standard error of the mean of individual results from 
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FIG. 2: (Color online) Fluorescence spectra (o) for 9,10,11 Be + 
in the 2s\/2 — * %Pi/2 transition as a function of the Doppler- 
tuned frequency in collinear (left) and anti-collinear (right) 
excitation. Frequencies are given relative to the respective 
hyperfine-structure center of gravity for the odd isotopes and 
the resonance frequency for 10 Be. Lines (— ) are fitting results 
for Voigt profiles as discussed in the text. 



all measurements on a particular isotope. An uncertainty 
of 0.5 MHz accounting for possible misalignments of the 
laser and ion beam overlap was added quadratically. Re- 
coil corrections in the isotope shift are only at the 10% 
level of the final uncertainty. Charge radii are deduced 
by subtracting the mass-dependent isotope shift of the 
respective isotope pair as calculated in Ref. 0. The re- 
maining nuclear-volume shift provides the change in the 
mean-square nuclear charge radius 6 (r%) between two 
isotopes. Absolute charge radii r c must be related to at 
least one isotope for which the absolute radius is known 
and can then be deduced according to 
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with the theoretically calculated electronic factor C = 
-16.912 MHz/fm 2 0]. 

For stable 9 Be, r c was determined by elastic elec- 
tron scattering (l7|, and by muonic atom spectroscopy 
0. Reported charge radii of 2.519(12) fm [T3] and 
2.39(17) fm differ by 0.13 fm but agree within the rather 
large uncertainty of the muonic atom result. We have 
used the electron scattering result, but we note that this 
was not obtained in a completely nuclear-model indepen- 



TABLE I: Isotope shifts Su^ in the Dl line and theoreti- 
cal mass shifts Sv^s 0> f° r A Be " 9 ^ e - Uncertainties for 
the absolute charge radius r c include the uncertainty in the 
reference radius r c ( 9 Be) = 2.519(12) fm [T3]. 
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FIG. 3: (Color online) Experimental charge radii of beryllium 
isotopes from isotope-shift measurements (•) compared with 
values from interaction cross-section measurements (o) and 
theoretical predictions: Greens-Function Monte-Carlo calcu- 
lations (+) Fermionic Molecular Dynamics (A) [HI, 
ab initio No-Core Shell Model (□) 0, 0, [Si . 



dent way. Hence, reanalysis of the Be scattering data, as 
performed for the proton 21 1, would be useful to deter- 
mine whether the small uncertainty of 0.012 fm is reli- 
able. However, a small change in the reference radius 
causes primarily only a parallel shift of all charge radii. 

Results for <5(r 2 ) and r c arc listed in Table [TJ Isotope 
shifts in the D2 line have also been measured and are 
still under evaluation. The extracted charge radii agree 
with the values reported here, but are less accurate due 
to unresolved hyperfme structures. 

The derived nuclear charge radii are shown in Fig. [3] 
r c decreases from 7 Be to 10 Be, but then increases for 
11 Be. The decrease is probably caused by the cluster- 
ization of 7 Be into an a- and a triton-cluster, whereas 
9,10 Be are considered to be a + a + n and a + a + n + n 
systems, respectively, and are more compact. Accord- 
ing to a simple frozen core two-body model the increase 
from 10 Be to 11 Be can be attributed to the motion of 
the 10 Be core relative to the center of mass. Using 
r;? ( n Be) = R 2 cm + r^( 10 Be) a rms distance of 7.0 fm be- 
tween the neutron and the center of mass of 11 Be can be 



extracted directly from S (V 2 )- 
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TABLE II: Magnetic dipole coupling constants A of the 2s\/2 
and the 2px/ 2 states and the magnetic moments of the odd 
isotopes. 



Isotope A?s 
Y Be -742.90(25) 



(MHz) 



A 2p (MHz) f4m) 



Rcf. 



'Be 



Be 



-624.97(4) 

-625.008837048(10) 

-2677.4(8) 



-140.17(18) 

-118.00(4) 
-118.6(36) 
-505.4(5) 



-1.3995(5) 
-1.398(15) 

-1.177432(3) 

-1.6813(5) 

-1.6816(8) 
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To test nuclear-structure theories, predictions from dif- 
ferent models are included in Fig. Reported point- 
proton radii r pp were converted to nuclear charge radii 
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r c by folding in proton R p [21| and neutron R n 
charge radii and adding the Darwin-Foldy term 
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Nuclear-model dependent charge radii extracted from 
interaction cross-section measurements [28| show a sim- 
ilar trend as our data but changes are smaller for the 
isotopes below 10 Be and overestimated for 11 Be. All 
theoretical model predictions show good agreement with 
our results. Depending on the model, specific features 
are especially well reproduced. In the case of GFMC 
calculations @, [22| the agreement for 9 ' 10 Be is striking 
whereas the decrease between 7 Be and 9 Be is slightly 
underestimated. Unfortunately, a charge radius calcula- 
tion for n Be has not been obtained to date. FMD [23} 
and improved NCSM [3, HI HI] calculations include the 
halo nucleus and the change between 10 Be and 11 Be is 
very well predicted. Concerning the other isotopes FMD 
does a better job for the isotope pair 9,10 Be, whereas the 
change in the pair 7,9 Be is best described by the NCSM 
results (note that the absolute charge radius depends on 
the reference radius). 

From fitting the observed hyperfine structures we ob- 
tain the magnetic dipole constant A for the 2s and the 
2pi/2 states as listed in Tab. Jill Based on the precision 
measurement of Ai s in 9 Be [30] and the corresponding 
nuclear magnetic moment \i = —1.177432(3) /zn [IH, we 
determine the magnetic moments of 7 Be and 11 Be. For 
11 Be we confirm the previously measured value from an 
optical pumping /3-NMR experiment |32l| . For 7 Be, an 
earlier value from optical hyperfine measurements [291 ] is 
improved by more than an order of magnitude. 

To summarize, we have measured the charge radii of 
7,9,io,iig e with on _li n e frequency-comb based collincar 
laser spectroscopy. From a simple frozen-core two-body 
model we obtain a rms distance of about 7 fm between the 
halo neutron and the center of mass in n Be. Comparison 
with elaborate nuclear structure calculations (GFMC, 
FMD, and LBSM) shows in all cases good agreement 



with our measurement. The magnetic moment of 11 Be 
has been confirmed and the accuracy for 7 Be consider- 
ably improved. The developed experimental technique 
has the potential to be applied to the very short-lived 
nuclei 12,14 Be and can generally improve the accuracy of 
collinear laser spectroscopy where this is required. 

This work was supported by the Hclmholtz Association 
(VH-NG-148), BMBF (06TU263I, 06MZ215, 06UL264I), 
and the EU (FP-6 EU PJI3-CT-2004-506065). M.K. 
was supported by the EU (MEIF-CT-2006-042114). We 
gratefully acknowledge technical support by M. Fischer 
and R. Holzwarth from Menlo Systems and by the 
ISOLDE technical group. 



[1] 
[2] 

[3. 
[4 

[5 

[6 

[-: 

[8 

[9 
[10 

[11 
[12 

[13 
[14 

[is; 

[16 
[17 

[is; 

[19 
[20 
[21 
[22 

[23 
[24 

[25 
[26 
[27 

[28 
[29 

[30 

[31 



I. Tanihata et al, Phys. Rev. Lett. 55, 2676 (1985). 

S. C. Pieper and R. B. Wiringa, Ann. Rev. Nucl. Part. 

Sci. 51, 53 (2001). 

R.Neugart et al, Phys. Rev. Lett. 101, 132502 (2008). 
Z.-C. Yan and G. W. F. Drake, Phys. Rev. A 61, 022504 
(2000). 

Z.-C. Yan and G. W. F. Drake, Phys. Rev. Lett. 91, 
113004 (2003). 

M. Puchalski, A. M. Moro, and K. Pachucki, Phys. Rev. 
Lett. 97, 133001 (2006). 

Z. C. Yan, W. Nortershauser, and G. W. F. Drake, Phys. 
Rev. Lett. 100, 243002 (2008). 

M. Puchalski and K. Pachucki, Phys. Rev. A 78, 052511 
(2008). 

G. Ewald et al, Phys. Rev. Lett. 93, 113002 (2004). 
L.-B. Wang et al, Phys. Rev. Lett. 93, 142501 (2004). 
R. Sanchez et al, Phys. Rev. Lett. 96, 033002 (2006). 
P. Mueller et al, Phys. Rev. Lett. 99, 252501 (2007). 
A. Csoto and K. Langanke, Nucl. Phys. A 636, 240 
(1998). 

P. Navratil, C. A. Bertulani, and E. Caurier, Phys. Rev. 
C 73, 065801 (2006). 

M. Zakova et al, Hyperf. Interact. 171, 189 (2006). 
T. Nakamura et al, Phys. Rev. A 74, 052503 (2006). 
J. A. Jansen, R. Th. Peerdeman and C. de Vries, Nucl. 
Phys. A 188, 337 (1972). 

R. Neugart, Nucl. Instr. Meth. 186, 165 (1981). 

W. Geithner et al, Hyperf. Interact., 129, 271 (2000). 

L. A. Schaller et al, Nucl. Phys. A 343, 333 (1980). 

I. Sick, Phys. Rev. Lett. B 576, 62 (2003). 

S.C. Pieper, K. Varga and R.B. Wiringa, Phys. Rev. C 

66, 044310 (2002). 

R. Torabi, GSI, priv. comm. (2008). 

C. Forssen, P. Navratil, W.E. Ormand, and E. Caurier, 
Phys. Rev. C 71, 044312 (2005). 

P. Navratil, priv. comm. (2008). 

S. Kopecky et al, Phys. Rev. C 56, 2229 (1997). 

J. L. Friar, J. Martorell, and D. W. L. Sprung, Phys. 

Rev. A 56, 4579 (1997). 

I. Tanihata et al, Phys. Lett B 206, 592 (1988). 

S. Kappertz et al. in ENAM98: Exotic Nuclei and Atomic 

Masses, B. M. Sherrill et al. (Eds.), AIP 110 (1999). 

D. J. Wineland, J. J. Bollinger, and W. M. Itano, Phys. 
Rev. Lett. 50, 628 (1983). 

W. M. Itano, Phys. Rev. B 27, 1906 (1983). 



[32] W. Geithner et al, Phys. Rev. Lett. 83, 3792 (1999). 



